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We present a fully relativistic analysis of Bessel beams revealing some 
noteworthy features that are not explicit in the standard description. It is 
shown that there is a reference frame in which the field takes a particu- 
larly simple form, the wave appearing to rotate in circles. The concepts of 
polarization and angular momentum for Bessel beams is also reanalyzed. 
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Bessel beams are solutions of Maxwell's equations with important proper- 
ties: they propagate with an intensity pattern that remains invariant along a 
given axis P, and carry an angular momentum that is not due to their polar- 
ization state j2] • Their experimental realization in laboratories has attracted 
much attention in recent years [3]. 

The aim of the present paper is to study the general properties of Bessel 
beams within a fully relativistic framework. We show that, contrary to the 
case of a plane wave, there is a privileged reference frame in which the linear 
propagation is eliminated and the wave propagates circularly. This particular 
frame corresponds, in a certain sense, to an antiparaxial limit. A general 
implication is that Maxwell's equations admit as solutions electromagnetic 
waves that propagate in circles. 

Our starting point is the definition of a Bessel electromagnetic wave in 
terms of Hertz potentials 111 and II2 jl]. In cylindrical coordinates {p, <f),z}, 
the electromagnetic potentials are given by [5] 

• = -|n„ (i) 

(in units with c = 1), and both Hertz potentials (i = 1,2) satisfy the 
equation: 

Any solution of this equation regular at the origin can be written as a linear 
combination of the functions 

n x = -r^Jm(k±p) exp{-iut + ik z z + zm0}, (4) 

K 1 



and 



n 2 = -r^J m (k±p) exp{-iut + ik z z + imcj)}, (5) 
k 1 



where J m is the Bessel function of order m, £ K and Bk are constants, and 



k± = ^Juj 2 — k 2 z is the transverse wave- number; here and in the following, 
the set of numbers {k^,m, k z } will be denoted with the generic symbol K. 
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An electromagnetic mode is associated to each Hertz potential, giving rise 
to transverse magnetic (TM) and electric (TE) modes respectively. Using a 
Cartesian basis, it follows that the most general superposition of these waves 
is given by [B] 



E 



K 



-iujt+ik z z 



2k i 



(u£ K + ik z B K )J m -i(k ± p)e 
+ (ujS k - ik z B K )Jm+i(k ± p)e 



i(m-l)4> 



i(m+l)(f>( 



G.7' ~~\~ ^67, 



+ e 



-iuit+ik z z 



B K Jm{k±p)e im ^e z 



(6) 



and 
B K 



-iujt+ik z z 



2fcj 



[ik z £ K - uB K )J m -i(k ± p)e 
{ik z £ K + uB K )J m+1 (k±p)e 



i(m+l)(j> , 



+ e 



-iut+ik z z 



(7) 



Clearly the electric an magnetic fields are interchanged under a duality trans- 
formation Bk — * £k and £r — * —&k- A compact way to write the above 
expressions is: 



E K ± iB K = (£ K T*B K )e- M * z 
uj ± k 



2k i 



+ 



2k 



(8) 



(here and in the following, the argument of all Bessel functions is fcj_p). 

The complex Lorentz invariant of the electromagnetic field follows from 
the above expressions: 

(E K + iB K ) ■ {E* K + iB* K ) = 



{£k — iBk){£k ~ iB*K)\\Jm+l\ 2 + |^m-l| 2 — 2| J n 



and also 



(jftE K ) 2 - (JRB K ) 2 + 2i{^E K ) ■ (m K ) 



(£ K - iB K )e 



-iLot+ik z z-{-im<f) 



+ (£* K - iB* K )e 



* \ „iojt—ik z z—imd 



(9) 
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+ -(S K - iB K )(£* K - iB* K ){J m+l - J m _!) 2 , (10) 

where 3?E and !RB are the real parts of the electric and magnetic fields, that 
is, the physically measurable quantities. Notice that, unlike the case of a 
plane wave, the Lorentz invariant 3?Ex • 3?B^ can be different from zero. 
From Maxwell equations, this implies that 

5KE X ■ V x UE K ^ (11) 

for monochromatic waves. A similar expression follows for the magnetic 
field. Eq. (II lj) can be satisfied only if there is no scalar field g(r) such that 
V x (g$tEi K ) = 0. Thus the rotational behavior of the electric and magnetic 
field cannot be removed continuously: this is a manifestation of the fact that 
Bessel beams have optical vortices. 

There are two particularly important superpositions of Bessel modes: the 
so-called right and left polarized states, given by the conditions 

u£ K ± ik z B K = 0, (12) 

and the states for which 

£ K ±iB K = 0. (13) 

The first condition corresponds to a wave whose electric field has components 
along either e x + ie y or — ie y , and its z component is negligible in the 
paraxial approximation; however, the corresponding magnetic field in this 
case has a rather complicated structure. The second definition appears more 
naturally within a quantum optical framework: it corresponds to a basis of 
orthonormal modes in which the operators for the energy, the helicity, and 
the z-components of linear and orbital angular momentum are all simulta- 
neously diagonal 0. Eqs. (fT2"j) or Eqs. (fE^) are equivalent in the paraxial 
approximation. 

For circularly polarized beams in the sense of Eq. (fT3*j) . the complex 
Lorentz invariant vanishes, in close analogy with a standard plane wave of 
arbitrary polarization, while the Lorentz invariant given by Eq. (jlLip takes a 
simple form but it is not strictly zero. This means that the temporal averages 
of |3?Ea'| 2 — |^Bft;| 2 and 3?E/< • KB a- are zero, although the corresponding 
values at a given time do not vanish. 

At this point, we notice that there is a particular reference frame in which 
TE and TM Bessel modes take a simpler form. This is the frame moving 
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along the z axis with velocity v = k z /u and Lorentz factor 7 = u/k±. The 
standard Lorentz transformation to this new frame is 



t' = - — (ojt — k z z) 



— (uz - k g t), 
k 



(14) 



and the same transformation holds for the t and z components of the potential 
A^. Thus the transformed electromagnetic field does not depend on the 
coordinate z', since k' z = 0. Changing to this frame is equivalent to taking 
an antipar axial limit. If the mode satisfies the paraxial approximation in 
the laboratory frame, the above Lorentz transformation involves an ultra 
relativistic velocity, v ~ 1. 

A direct calculation shows that the scalar potential $ = in the new 
frame, and therefore the Coulomb gauge is satisfied directly. As for the 
vector potential, it takes the form: 



A' = - — {SkVLk + Bkv k ), 
k± 



(15) 



where 
and 



m 



k±p 



-ik±t'+im(j) 



Jm—1 ^ 



i(m— 1) 



(e x + ie y ) + J m+ i e 



267, 



-ik±t' +im<p 



e 2 . 



The electric and magnetic fields are 



E' = S K u K + B K v K 



and 



B' 



-B K u K + £ K v K . 



-^*'(16) 
(17) 

(18) 
(19) 



In the antiparaxial reference frame, the magnetic (electric) field of a trans- 
verse electric (magnetic) mode is parallel to the z axis. 

The standard definition, Eq. ([12JL of right and left polarization imple- 
mented in this frame leads to E-k = because k' — 0. That is, the "circularly 
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polarized beam" becomes a TM mode whose electric field points in the z di- 
rection. On the other hand the definition Eq. (J13|) leads to an electric field 
with a projection in the XY plane that is a superposition of two opposite cir- 
cular vectors with amplitudes proportional to J m=F i, while the z component 
of the field is proportional to J m and has a relative phase ±7r/2. 

In Fig. (1), we illustrate the electric and magnetic fields in the plane 
perpendicular to the z axis, as seen in the antiparaxial frame, for a circu- 
larly polarized mode in the sense of Eq. (|13|). Notice that, as previously 
mentioned, 3?E and 3?B are not perpendicular in general. Edge phase dislo- 
cations, saddle points and vortices [7] are present. Some of these structures 
are formed around the zeros of J m and J' m . In Fig. (2), the intensity patterns 
1 3F£E | 2 , |3?B| 2 and their sum are illustrated. Notice the complementary space 
distribution of the electric and magnetic fields. 

As for the dynamical properties of the field, a noteworthy feature in the 
new frame is that the z component of the Poynting vector is 

(ME' K x m' K ) z = 1 -{S K B* K - B K Z* K ){Jl_ x - J 2 m+l ), (20) 

so that it vanishes for pure TE or TM modes, as expected from the fact that 
k' z = 0. On the other hand, the local flux of energy for circularly polarized 
modes in the sense of Eq. f)13|) is different from zero even in this frame. 
Nevertheless, the integral over whole space of {3tE' K x ^RB' K ) Z is zero since 

poo 4771 f°° 

/ [J m +i{k±p) - J^ l _ 1 (k±p)]pdp = — J' m (k±p)J m (k±p)dp = 0. (21) 
Jo k± Jo 

Let us now turn our attention to the angular momentum. Its total density 
is defined as 

^=^ rX ( ExB )' ( 22 ) 

and it is known that in the Coulomb gauge, up to a surface term, it can be 
decomposed into the sum of the so-called orbital angular momentum density 

£ = ^E E *( rx V)A, (23) 

i 

and the spin angular momentum density 

S = — E x A 

47T 
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(24) 



The above equations can be obtained in a relativistically covariant form. 
Since the conservation of angular momentum is related to a rotational sym- 
metry, it is natural to use the Killing vector k a associated to rotations around 
the z axis. This vector is defined as k a d a = d/d<fr, and has the standard 
property V ' a kp = —Vpk a) where V a is the covariant derivative (see, e. g., 
Weinberg jU]). In cartesian coordinates = to 3): 

k a = (0,-y,x,0) , V a ^ = e to/J3 • (25) 
Using the energy- momentum tensor for the electromagnetic field F^: 

T a P = — \F a »F p ■ - -g a P(F x >*F Xtl )} , (26) 



An 
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it follows that the condition VpT alS = implies that the four-vector J a = 
kpT 01 ^ is conserved, that is: V a J° = 0. Thus, using the fact that F a p = 
VaAp — V ' pA a and V ' a F a/3 = 0, it turns out that J a can be written as: 



47T 



F^kTVvAn + F^A^V^K - V^F a »k v A v ) - ^k a (F x »F x ^)} . (27) 



Since J a has zero divergence, the integral over a three-dimensional hypersur- 
face with normal unit four- vector n a and volume element dV, 

J = J dVn a J a , (28) 

is independent of the particular choice of the hypersurface and thus is a con- 
served quantity. In particular, choosing such hypersurface as t = constant, 
which implies that n a = (1,0), it follows that 

n a J a = - 1 - \F 0m k»V u A m + F^A v V^k u - V m (F° m k u A u )} . (29) 

We readily identify the first term in this equation as the orbital angular 
momentum density, given by Eq. (|2H|) and the second term as the spin den- 
sity, Eq. (J23J). As for the third term, its volume integral involves a three- 
dimensional divergence and can be taken as zero if the field vanishes at 
infinity. Thus, we have recovered the standard formulas ()23|) and ()24j) in a 
covariant form. The dependence of the reference frame appears through the 
choice of time-like unit vector n n . 
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Applying Eq. (J23)) and taking only the real parts of the complex fields, it 
turns out that in the laboratory frame 
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K 



in 
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u, 2 \£ 



K\ 



K\B 
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m— 1 



J, 



m+1 



+ 2 u 2 \£ K \ 2 cos2(ip + ip E ) + k 2 z \B K \ 2 cos 2(<f + ip B ) J m _iJ m+ i 
+ 2k z u\£ K \\B K \sm((p B -ip e )(J i 
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I 2 N 



2k 2 ± \B K \ 2 [l + cos2( V + VB )}Jl} 



(30) 



where tp = mcf) + k z z — ut, and tpg and ps are the phases of the complex 
amplitudes £k and Bk- The density averaged over a cycle is accordingly 
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K 



m 



167TLdk 2 t 



u 2 \£ 



K\ 



k 2 \B 



K\ 



r 4- T 



m+1 



+ 2k z u\£ K \\B K \ sin(^ B - ipe)(J^-i - J^+x) 
+ 2k\\B K ~ 



1 2 t2 



(31) 



which is non null for all kinds of polarization and is proportional to the az- 
imuthal number m as expected. As for the spin density of the electromagnetic 
field, it follows from the previous equations that 



•■it 



IQiiujk 2 ^ 
+ 2k z u\£K 



u 2 \£ 



K\ 



k 2 \B 



K\ 



m+1 



7 2 N 



\B K \ sin(v9 jB - Pe){J\ 



m+1 



I 2 ' 



(32) 



Thus, in the antiparaxial frame a purely transverse magnetic mode has null 
Sf, while any other cylindrical mode, such as a circularly polarized wave, 
has a "spin" density in the z-direction. Notice, however, that the integral 
of the z component of the Poynting vector, Eq. (J2(Jj) . and the integral of Sf 
over the whole space are zero in this frame. This result is consistent with the 
interpretation of the integral of Sf as a helicity operator. As for the sum of 
the orbital and spin angular momentum densities averaged over a cycle, it 
turns out to be 



A 



\<OT\k 



{\£ K \ 2 [(m - l)Jl_, + (m + l)J 2 m+1 ] + \B K \ 2 mJ 2 m } (33) 



in the antiparaxial frame. 

An alternative definition of the orbital and spin angular momentum den- 
sities that has the advantage of being gauge independent has been proposed 
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by Barnett^Ul, who defined the time-averaged z-component of the orbital 
angular momentum flux (for complex electric E and magnetic B) as 



M. 



ORB 



and the spin flux as 



b: 



dB* 



dB* 



M. 



SPIN 
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iTXUJ 



)> — i{E x B* + E y B*) 



Applying these formulas to a Bessel mode we find that 

k 2 



M. 



ORB 

zz 



and 



+ 



2{k z /kj 



\£k\\Bk\ • 



B: 



dE 3 



(34) 



(35) 



{\£ K \ 2 + \B K \ 2 ) [(m - 1) J^_! + (m + 1) J, 



l^k\ 
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16tiuj 

sin(v9 £; - ^ B )[(m - 1) J^_ x - (m + 1) J m+1 ] 



2 

m+1 



(36) 



SPIN 
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I^A'I )(<^ m _l J m+1, 



2{k z /k A 



\£k\\Bk\ 



m+l t 



(37) 



Thus, according to this definition, the total density of angular momentum 
averaged over a cycle is 



M, 



Mg RB + M. 
mk, 

£ 



SPIN 



m 
IQttuj 



K 



B K \ ){J m ~i + J, 



2 

m+l 



2{k z /k ± ) 2 + 1 \S K \\B K \ sin(^ - <f B ){J m ^ - J m +i)(38) 



As a consequence, in the antiparaxial frame, both M® Z RB and M^f IN van- 
ish for pure TE or TM modes, but not for circular modes. This appar- 
ent inconsistency may be due to the fact that the total M zz transforms as 
M zz — > M zz +-/vJ z under a Lorentz boost along the z axis, unlike J z which is 
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invariant, so that its value depends on a given reference frame moving along 
the beam. 

In conclusion, a relativistic analysis reveals the existence of a particu- 
lar moving frame in which Bessel beams have a simpler form. In practice, 
this may be a frame moving at relativistic speed. However, since we are 
dealing with exact solutions of Maxwell's equations, some invariants can be 
easily calculated and the covariance of variables such as the angular mo- 
mentum densities can be studied. Although the definitions of the so called 
orbital and angular momentum densities have some ambiguities, it must be 
expected that these quantities do not change qualitatively, or even vanish, by 
a Lorentz transformation. Actually, we have shown that there are important 
cancellations in the antiparaxial limit. Thus it seems more appropriate to 
interpret the standard angular momentum density Sf as the helicity den- 
sity, in accordance with our discussions in a previous paper [3]. Moreover, 
the orbital angular momentum is expected to be finite in the antiparaxial 
frame if it is different to zero in any other frame, but that is not the case 
when the expression of M® Z RB is applied to TE or TM modes. It is also 
clear from our analysis that the definition of polarization which is used for 
plane waves cannot be applied unambiguously to a Bessel beam. One must 
be careful in defining what is meant by polarized states: in fact, we have 
shown two different definitions leading to different properties. As a further 
application of the present relativistic analysis, we will study the motion of 
charged particles in the field of Bessel beams in a forthcoming paper. 
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1 Figures captions 

Figure 1. Magnetic and electric fields in the plane perpendicular to the z axis, 
as seen in the antiparaxial frame, for a circularly polarized mode [defined by 
Eq. jnj)] with m = 2. 

Figure 2. Intensity patterns |3?E| 2 , |3?B| 2 and their sum, as seen in the 
antiparaxial frame, for a circularly polarized mode [defined by Eq. 1)130] with 
m = 2. 
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